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A THEOREM ON CENTRAL VELOCITY DISPERSIONS 

Jin H. An,''^ and N. Wyn Evans' 
j inUdark-cosmology . dk , nweOast . cam . ac . uk 

ABSTRACT 

It is shown that, if the tracer population is supported by a spherical dark halo with a core or a cusp diverging 
more slowly than that of a singular isothermal sphere, the logarithmic cusp slope y of the tracers must be given 
exactly by y = 2y6 where /3 is their velocity anisotropy parameter at the center unless the same tracers are 
dynamically cold at the center If the halo cusp diverges faster than that of the singular isothermal sphere, the 
velocity dispersion of the tracers must diverge at the center too. In particular, if the logarithmic halo cusp slope 
is larger than two, the diverging velocity dispersion also traces the behavior of the potential. The implication 
of our theorem on projected quantities is also discussed. We argue that our theorem should be understood as a 
warning against interpreting results based on simplifying assumptions such as isotropy and spherical symmetry. 
Subject headings: stellar dynamics — galaxies: halos — galaxies: kinematics and dynamics — dark matter 
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1. INTRODUCTION 

Gravitationally-interacting coUisionless N-body systems 
are a good model for both dark-matter halos and stellar sys- 
tems, and therefore have been a subject of many studies. Early 
on, it was realized that the evolution of such systems is gov- 
erned by what is now known as the col lisionless Boltzmann 
equation (lEddingtonlll9Tl iJeansllloTl) . However, smce It 
deals with the distribution function, which is generally inac- 
cessible in all but very extensive and nearly complete data 
sets, its utility is somewhat limited in reality. An alterna- 
tive to dealing directly with the distribution function is to fo- 
cus on the relations among the statistical moments of it. The 
set of their governing equations is obtained by taking veloc- 
ity mome nt integrals on the collisionless Boltzmann equation 
(see e.g.. lBinnev & Tremainei2 008) and bears the name of Sir 
James H. Jeans (1877-1946). Of a particular interest among 
the set is the one resulting from the first moment integral, 
which is what is usually referred to as the Jeans equation. The 
equation summarizes the momentum conservation in a local 
volume element and is an analogue to the fluid/gas-dynamical 
Euler equation except for the presence of anisotropic stress 
tensor relating to the local velocity dispersions in place of the 
usual pressure term. 

In recent years, interest in the Jeans equations has 
been rekindled for various reasons. First, N-body simu- 
lations exhibit striking regularities among dynamically re- 
laxed dark matter struc tures (e.g., iTaylor & Navarrol 1200 II: 
[Hansen & Moord l2006h . It is reasonable to suspect that 
these are the results of some underlying physics working 
to erase the memory of the initial conditions and settle the 
structure into a 'universal' form constrained by the Jeans 
equat ions (e.g., iDehnen & McLaughUnI 120051: iBarnes et"aLl 
|2006|) . Second, a wealth of new data has become available 
on the velocities of giant stars in neM"by dw arf spheroidals 
dKlevna et all 120021: IWkinson et al.l |2004 and of plane- 
tary nebulae in nearby e lhptical galaxies lRomanowskv et al.1 
I2003t iDekel et al.ir2005h . The structure of dark halos can be 
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mapped out through these data coupled with the Jeans equa- 
tions. Of course, it is of great interest to establish the dark 
halo structure in the central regions, whi ch is predicted to be 
cusped in hierarch ical cosmologies re.g.. lNavarro et al.lll995t 
lMooreetal.llI998l) . 

All these have motivated greater theor e tical s crutiny of the 
Jeans equations. Recently. lEvans et al.l (l2009l) have shown 
that any cusped profile of an isothermal tracer population can 
be supported by the potential generated by a dark halo only 
if the halo possesses an isothermal cusp. In this paper, we 
extend their study and subject the properties of systems gov- 
erned by the spherical Jeans equations to a thorough theoreti- 
cal analysis. 

2. JEANS EQUATIONS IN SPHERICAL SYMMETRY 

The Jeans equations under the spherical-symmetry and the 
steady-state assumptions reduce to 



■2/3- 



vcr. 



d^ 



= -V- 



(1) 



dr r dr 

where v = v(r) and cr,. = cr,.(r) are the density profile and 
the radial velocity dispersion of the tracer population. The 
density here is nominally the number density, which follows 
the derivation of the Jeans equations from the collisionless 
Boltzmann equation. However, since equation ([T]l is linear 
in V, it is still valid for any constant multiple of v, such as 
the mass or the luminosity density, provided that the tracer 
population has homogeneous properties. 

The velocity anisotropy parameter /3 = /3(r) is defined such 
that 



(2) 



where o-g{r) is the 1-d tangential veloc ity dispersion of the 
same tracers dBinnev & Tremain3l2008l) . The luminous trac- 
ers are moving in a gravitational potential ^(r), which, 
through the spherical Poisson equation. 



r^ dry dr j ' 



(3) 



is related to the density profile p(r) of the dark halo. Strictly 
speaking, p in equation O) is the total mass density that in- 
cludes all gravitating masses. Here, it will be simply re- 
ferred to be a dark halo, which is basically a label that sig- 
nifies that we do not demand that the potential should be self- 
consistently generated by the density profile of the tracers. 
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Although physical solutions are still subject to constraints that 
they must be non-negative and the tracer mass density may not 
be greater than that of 'dark halo' anywhere, these will not be 
considered explicitly in this paper. 

3. POWER-LAW SOLUTIONS TO JEANS EQUATIONS 

We begin our study with an analysis of power-law solu- 
tions to the Jeans equations with constant anisotropy param- 
eter. These are often good approximations locally, and have 
the advantage that they are analytically tractable. 

The spherical Jeans equation ^ is always formally inte- 
grable such that 



where Q = Q{r) is the integrating factor 

InQ = Z'^'"^- 
If V and 'I' behave locally like a power law. 



(4) 



(5) 



dr 



B 



where A,B > Q and 6 < 1, we can find solutions to the Jeans 
equations once the behavior for the anisotropy is prescribed. 
If the potential is self-consistently generated, the power in- 
dices and normalization constants are related to each other 
such that 7 = 5 H- 2 and AnGA = B{\ - 6) (see eq. ID). 

The easiest assumption to make regarding the behavior of /? 
is that it is constant. In reality, this probably is not true, but 
here we are interested in the generic behavior of solutions for 
given local power-law assumptions on v and Therefore, as- 
suming constant /3 is valid provided that its variation is much 
slower than that o f the density and the potential. Notably, 
[Hansen & Moord (l2006h found that the logarithmic density 
slope and the anisotropy parameter are linearly related in sim- 
ulated dark halos. Although thek detailed findin g appears not 
to be always true (see e.g., iNavarro et al.ll2008h . the general 
idea seems to be still valid. That is to say, the spatial variation 
of the anisotropy parameter is sufficiently slow so that it can 
be considered to be locally constant while the density profile 
is approximated as a power law. The situation is believed not 
to be much different in stellar systems, for which no evidence 
to the contrary is obvious. 

Under the assumption of the constancy of j3, we have Q - 
r^^, which results in 



B 



X+Cr 



y-213 + 6 
r^-^P {B\nr-^ + C) 



if 6 2/3- r, 
if5 = 2j8-7. 



(6) 



where C is an integration constant to be determined from the 
boundary condition. 

While the solution in equation (|6]l is always valid for given 
assumptions, it is somewhat easier to follow its behavior if we 
consider a finitely-deep and an infinitely-deep potential well 
separately. In the next two subsections, we investigate each 
case in detail, and find the constraints on the behavior of the 
radial velocity dispersion provided by the Jeans equations. 

3. L Case 1: finite central potential wells 

If the dark halo diverges like a singular isothermal sphere 
(which behaves as r"^; henceforth SIS), the resulting potential 



is logarithmically-divergent and so 5 = 0. Thus, given the 
power-law assumption, for any cusped halo diverging slower 
than a SIS, we can limit 6 < 0. Letting p - -5 > Q (that is, 
the potential behaves like ^ "Fq -t- Br''/ p), the leading term 
for crj as r —> becomes 



B 



rP{B\n 



\fp<y- 2/3; 



r-i + C) 



ifp 
if y ■ 



■7-2/3; 
2/3 < p. 



(7) 



For < p < y - 2/3, we have cr^ ~ r'' —» as r —> 0. If 
p - y - 2p > 0, then cr^ ~ In r"' and again limr^o o"? — > 0. 
Moreover, the logarithmic slope of cr^ still tends Xo p > Q 
in the limit of r — > 0. Finally, \f y - 2/3 < p, then cr^ ~ 
f-y-^p^ However, for this last case, if y < 2/3, then — > oo 
as r — » 0. Since any finitely-deep central potential well is 
unable to support tracer populations with divergent velocity 
dispersions, the physical possibilities are limited to be y > 2/3 
(c.f. An & Evans 2006). 

In conclusion, with a finite central potential well, the possi- 
bilities are either (i) cr^ — » as r — > [here, the logarithmic 
slope of crj tends to mm(p,y - 2/3) > 0] or (ii) y - 2/3 with 
finite and non-zero cr^ at r = (see also lEvans et alJl2b09h . 

3.2. Case 2: centrally divergent potentials 

Next, we consider a centrally-diverging potential, for which 
Q < 6 < I. Here, the 6 - Q case corresponds to a logarithmic 
potential and a SIS-like dark halo cusp, whereas a point mass 
potential is represented hy 6 - I. 

For these cases, the theorem of lAn&EvaiisI (l2006h pro- 
vides us with the constraint that y > ^6 + /3(2 - 6), and so 
that 6 -2/3 + y > (| - /3)6. Since /3 < 1, we finally find 
that 5 - 2/i + y > 0. Here, this is strictly larger than zero if 
6 0. Consequently, the leading term for cr^ as r — > with a 
divergent potential is 



B 



+ Cr 



■7-2/3 



y-2/3 
Blnr' +C 

B 



1 



if 5 = and r > 2/3; 
if 6 ^y -2/3^0, 



y-2/3 + 5 r^ y-2/3 + 5 



(8) 



(9) 



if < ^ < 1, for which y - 2/3 + 6 > 0. Here, note that 
Blnr-Mf 5 = and -{B/6)r-^ for < 5 < L In 
addition, p oc r^(2+s) ^^j. Q < 5 < 1 if the potential is generated 
by p, whereas B = GM, and 6 - I if a black hole of mass M, 
dominates the potential. 

In conclusion, cr^ in a divergent potential well traces the po- 
tential except for the case of a logarithmically-divergent po- 
tential with y > 2/3, for which it is non-zero and finite. 

4. PROPERTIES OF SOLUTIONS 
TO JEANS EQUATIONS AT CENTER 

The preceding results are interesting, but the arguments 
leading to them are restricted to power-law solutions of the 
Jeans equations with constant anisotropy parameter. In fact, 
the results hold good more generally albeit in a weaker form. 
In the following, we shall derive a general theorem that ex- 
tends the preceding results. In this section, we will set up the 
framework and proceed to the case concerning the halo with a 
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core or a cusp that is shallower than the SIS, which is likely to 
encompass most astrophysically interesting models. We will 
extend the theorem to all physically allowed models of halos 
and resulting potentials including those dominated by a cen- 
tral point mass in § |5] The theorem is deduced by analyzing 
the Jeans equation dU in the limit of r — > without refer- 
ence to constancy of the anisotropy parameter or power-law 
behaviors. The resulting constraints however are only strictly 
applicable to the central limiting values. 

4. 1 . Preliminaries 

First, let us recast the Jeans equation into a more useful 
form. We begin by integrating equation (O, which leads us to 



r^— = GM(r) = GM. + 4nG 



r 2 
I drr p(r) 

Jo 



... . (10) 

dr Jo 

where M(r) is the enclosed mass within the radius of r, and the 
integration constant M, represents a central point mass (e.g., 
a supermassive black hole). However, we postpone detailed 
consideration of the M, + case until §|5] 
Equation ([T]) is then equivalent, with the enclosed mass, to; 

GM{r) 



(11) 



Also we have introduced the logarithmic slopes of the tracer 
density and the velocity dispersion, namely (note the signs) 



r = - 



d\nv 
d\nr 



r dv 
V dr 



dlno-j _ r do-j ^^^^ 
dlnr crl dr 



In the following, we consider the behavior of the system at 
the center, as indicated by the limit of equation (fTTT i as r — » 0. 
All the subsequent arguments operate under the assumption 
that every quantity considered here is well-behaved, continu- 
ous and smooth. 

4.2. Systems with vanishing M/r at the center 

Here, we ba sically repeat the argument found in the sec- 
tion 5 of Evans et alj (|2009|) with a slight refinement. The 
result will form a part of the theorem to be proven in § |5] and 
highlights the astrophysically relevant information. 

The condition for the left-hand side of equation ( fTTT i to van- 
ish in the limit r — > is given by limr^o M(r) = M, = 
and, from I'Hopital's rule, dM/dr\,~Q - 0. The last bit is 
equivalent to lim,-^o pr^ = - that is to say, limr^o p is fi- 
nite (i.e., a cored profile) or p diverges at the center slower 
than a SIS. Hence, assuming M, = (i.e., no central point 
mass), if limr^o pf^ - 0, then the right-hand side of equa- 
tion (fTTT) should also vanish as r — > 0. This is possible only 
if (i) cr^Q = 0, or (ii) - jq - 2/3o- Here and throughout, 
the subscript "0" is used to indicate the limiting value at the 
center. 

Suppose that (ii) is the case. Here, if yo < 2/3q, then it would 
be that crj ~ r^'^A-ro) _> oo as r -> 0. However, the SIS, for 
which lim,-^o pr^ is non-zero and finite, can only generate a 
potential diverging as fast as logarithmic. Thus, the velocity 
dispersion that diverges like a power law cannot be supported 
by the potential generated by the density cusp that is shallower 
than that of the SIS. Consequently, the corresponding case, 
70 < 2y6(), is unphysical and not allowed. On the other hand. 



if 70 > 2y6o, then crj ~ r' 



.(yo-2/So) 



0, which reduces to the 



case (i). In conclusion, a spherical dark halo with a core or a 
milder cusp than that of a SIS (i.e., lim^^o pf^ - 0), can only 
permit tracer populations satisfying the constraint ySo - 7o/2 
or those with ct^q = 0. 



4.3. Vanishing central velocity dispersions 
The preceding discussion indicates that ct^q = is neces- 
sary for lim^^o pr^ = if 70 2y6o- The initial impression of 
equation ( fTTT i notwithstanding, ct^q = alone however is not 
sufficient for vanishing M/r, either Formally, this is because 
the behavior of equation (fTTT i cannot be specified for tracers 
for which ct^q = and ySo - -c^ without reference to the 
speed of each approach to its limiting value. 

From its definition (eq. ||2l), p = -00 if cr^ =0 unless 
cTg = too. Here, j6o = indicates a tracer population with 
purely circular orbits toward the center Note that in princi- 
ple it is always possible to construct any spherical model with 
purely circular orbits (alth ough such models ar e subject to res- 
onant over-stabilities; see lPalmeretal.lfT989h . The non-zero 
tangential velocity dispersion here is the result of the random 
orientations of the orbital planes while the circular speed is 
uniquely specified by the enclosed mass. That is to say, the 
local dispersion of the speed of the tracers is actually zero 
although the tangential velocity dispersion may be non-zero. 

On the other hand, for all finite values of fio, then matters 
are simpler. Since 



<o = (l-^o)<o = 0, 



we easily find that cr^y = indicates that the total 3-d velocity 
dispersion at the center also vanishes. That is, the system must 
be dynamically cold at the center 

The conclusion of § l4.2l therefore may be rephrased as fol- 
lows: in a spherical potential well generated by a halo that is 
cored or cusped less severe than a SIS, the only allowed popu- 
lations of tracers are those either consisting of purely circular 
orbits toward the center or exhibiting vanishing central veloc- 
ity dispersions unless the limiting values of the cusp power 
index 70 and the anisotropy parameter y6o at the center are con- 
strained such that 70 - 2/3o- 

4.4. Discussion 

At first glance, the result of § 14.21 may appear to be coun- 
terintuitive as though it seems to suggest no pressure support 
at the center whereas the Jeans equation is supposed to bal- 
ance the force. This reasoning is faulty because the actual 
'pressure' in this case is given by pcr^, not cr^. Given the den- 
sity cusp, the system possesses non-vanishing kinematic pres- 
sure at the center even though it is dynamically cold. This is 
obvious in the power-law solutions to the Jeans equation of 
§|3l which show that vcr'j ~ ,^™(p-y-W)^ For a system that 
is radially-biased or is otropic toward the center (i.e., /3 > 0, 
for which 7 > from lAn & Evansll2006l and so tracers with 
a hole at the center are not allowed), the central 'pressure' is 
therefore always non-vanishing as r — > (it would be actu- 
ally divergent unless = < p — y). The 'pressure' of the 
tangentially-biased system (jS < 0) on the other hand can be 
vanishing if p > j. However, in this last scenario, it can be un- 
derstood that the tangentially-biased system is preferentially 
rotationally-supported toward the center 

If cr^ indeed vanishes at the center, this implies that there 
are no radial orbits in the model and that the distribution 
function has the property lim/^2^() /(L^, £) — > 0. This seems 
unusual, but there are mechanisms known that can depopu- 
late the radial orbits - for exa mple, scattering by the central 
cusp (Gerhard & Binnev 1985) or the radial orbit instability 
tPalmer & Papaloizou 1987) - albeit at the cost of driving the 
system away from sphericity. 
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The result can also be applied to a self-consistent system (or 
equivalently interpreted as a constraint on the central velocity 
dispersion of the dark halo itself). In a spherical dark halo 
that has a core or a cusp less severe than a SIS, the central 
limiting value of the anisotropy parameter must be exactly 
half of the numerical value of the cusp slope unless the central 
velocity dispersion vanishes. In particular, a cored halo must 
have an isotropic velocity dispersion at the center if it is not 
dynamically cold there. 

5. THE GENERAL THEOREM 

So far, we have investigated some astrophysically important 
subcases. Here, we derive and prove the general theorem that 
makes precise the interlocking constraints between the central 
limiting values of the density and velocity dispersion of the 
tracers and the potential. Those who are primarily interested 
in the result should skip to § 15.21 The next subsection provides 
a rigorous mathematical analysis of equation ( fTTT i that leads to 
our result. 

5.1. A derivation of the theorem 

We start by giving the binary relation '~' its precise math- 
ematical meaning. In the following, it is understood to be the 
short-hand notation such that a ~ as r — > if and only if 
both limr^o(fl/^) and limr^o(^/fl) are finite. We also define 
the binary relations > and similarly. That is to say, a > b 
(or a < /?) as r — > if and only if lim,.^o a/^ is divergent (or 
zero). In addition, it is also to be understood that the limit is 
taken to be always r — > unless the explicit reference to the 
limit is given to override. 

Next, we consider the behavior of the left-hand side of 
equation (fTTT i in relation to the mass density profile that gen- 
erates the potential. From the argument of § 14.21 we find that 
Mjr decays to zero in the limit r — > if M, - and p < r^^. 
The corresponding potential is either finite at the center or di- 
verges strictly slower than logarithmic (i.e., 'I' < Inr). On the 
other hand, it attains a non-zero finite limiting value at r — > 
if and only if M{r) ~ r. This is equivalent to p ~ r^^ and also 
'I' ~ Inr. Finally, M/r diverges as r — > if M, is non-zero 
(implying the presence of a central point mass) or p > r^^ 
in the same limit. The corresponding potential also diverges 
but strictly faster than logarithmic (i.e., Y > Inr) in the same 
limit. 

The behaviors of M/r and *!* in relation to each other and 
their respective logarithmic slopes may be explored in further 
detail. Let us first define p, the logarithmic slope of M/r, i.e.. 



P = 



dMM/r) 
dlnr 



1 d (GM\ 
d'i'/drd'r[~) 



(13) 



Note that in the limit of r — > 0, the logarithmic slope of the 
potential (difference) necessarily tends to the same value as 
P(). In particular, if ^^(O) = is finite, it naturally follows 
that GM/r = r(d^/dr) and so 



lim 



t/lnl^y-fol 
dlnr 



lim 



GM/r 



■ PQ > 



using I'Hopital's rule. If limr^o ^ 
we also find that 



^>-^'o 

■oo on the other hand. 



r ^'"I'^'l r GM/r 
lim lim — — — = po ^ 0. 



^0 dlnr 



^0 ^> 



Although the usual conditions for FHopital's rule for this case 
are only strictly met if M/r diverges, THopital's rule can in 



fact be proven only assuming a divergent denominator. Thus, 
the result holds even though M/ r tends to a finite value (in- 
cluding zero) as r — > 0. Moreover, it is obvious that po = 
if — > oo and M/r is finite as r — » 0. Next, if po + 0, 
then M/r ~ \A^\ where A*!* = T for limr^o *P - or 
A*? = I*!* - for Y(0) = being finite. By contrast, 
that po = indicates that M/r lA*!*]. If the central poten- 
tial is additionally finite (e.g, T - To ~ [Inr"']"', for which 
M/r ~ [Inr]"^), then M/r 0. The behavior of M/r as 
r — » for a divergent central potential depend on how fast 
diverges relative to logarithmic divergence (In r"') - e.g., for 
1^1 ~ In Inr-', m ~ Inr, and |T| ~ j[lnrf, we have that 
M/r ~ [Inr"']-', M ~ r and M/r ~ Inr"', respectively. 

Next, we proceed to analyzing the right-hand side of equa- 
tion (fTTI) . Here, we do not consider the /3 - -oo case, which 
represents the formal possibility of building the system with 
purely circular orbits. Then, since M and crj must be non- 
negative, yo - 2/3() - ffo > 0. Moreover, yo ^ 3 from the 
constraint that the central mass concentration must be finite. 
Together with the assumption that jSq is finite, we find that 
(yo - 2j6o - Co) is divergent only if ao diverges to negative 
infinity. However, then cr^ diverges faster than any power law 
to an essential singularity (e.g., e'^''). This is physically im- 
possible, because no real potential diverges faster than 1 /r nor 
is thus able to support such steeply diverging velocity disper- 
sions. Therefore, we limit yo - 2ySo - oro to be finite. 

If yo - 2j6o - (To > 0, it is clear that cr^ ~ GM/r, which also 
indicates that ao - po. If oq - po + 0, then cr^ ~ lA^F], too. If 
tto = = on the other hand, the behavior of (j\ still traces 
that of M/r, but cr]. ~ M/r < \A^\. 

If yo - 2jSo - ao, then (t\ > M/r and so ffo < po. Here, if 
the central potential is finite, then we find that < gp = yo - 
2j6o < po, from the constraint of I An & EvansI (l2006h . For a 
divergent potential (introducing 5o = -po, for which < 5o ^ 
1). the constraint of lAn & EvansI ( |2006|) . m > i(5o+A)(2-(5{)), 
indicates that \6q - /3oSo < yo - 2y6o = ffo < po - So- 
Now, if ^0 > 0, this would imply ySo > |. This is obviously 
impossible, and therefore 5o = po - ao - yo - 2j6o - 0. In 
addition, if M > r, it is clear that cr^ — > oo. Furthermore, from 
equation ( fTTT) recast to be 

r— = (y-2yS)cr -r— ^, 
ar ar 

we find for sufficiently-fast-decaying j - 2/3 that cr^ ~ \'¥\. 
This essentially implies that cr^ cannot diverge faster than 

5.2. The statement of the theorem 

In summary, the spherical Jeans equations permit only re- 
stricted physical possibilities regarding the limiting behaviors 
at the center. In particular, the central limiting value of the 
velocity anisotropy (/So; eq. 121) and those of the logarithmic 
slopes of the luminous tracer density (yo; eq. lfT2l ). the ra- 
dial velocity dispersion (ao; eq. lfT2l ) and the potential (po; 
eq. |[T3l ) must meet one, and only one, of the following list of 
choices, 

(i) pa = Q-o < yo - 2y6o and cr? ~ M/r, 

(ii) Po ^ ao - yo - 2y6o > and M^o is finite, 

(iii) Po = ao = yo - 2y6o = and lim,.^o *P = 

(iv) ySo = -oo. 
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Focusing on the behavior of the velocity dispersion, the resuh 
with the proviso fi > -oo is summarized as 



Hm 



d In crj. 
dlnr 



min(2 -Tcyo ■ 
-(To - 2) < ro - 



2ySo) > 
2A) 



(r,) < 2) 

(r,) > 2) 



(14) 



although this does not include all the information encom- 
passed in the above choices. Here, 

^ ^ d\np 
dlnr 

so that p - 2 - F, and extending to include the central point 
mass by setting Fq = 3. 

For a prescribed behavior of M/r or ^F, the above list re- 
turns the natural extension and generalization of our earlier 
results. If M/r — > for example, then either I) crj — > with 
crj ~ M/rorao = jQ-lfio > 0, or 2)0-0 - To- = 0- Con- 



sequently, we recover the conclusion of § 14.21 The implication 
of the list however is more detailed. First, if jo-lfio > po > 
(po > is necessary for vanishing M/r), then only the case (i) 
is possible and so cr^ ~ M/r and q-q - pQ. If po > addition- 
ally, then crj ~ - >Pol with a finite whereas crj < {A'i'l 
for Po = 0. On the other hand, with po > Jo - 2y6o > 0, we 
have ao - Jo - 2/3o > and so cr^ —» (and po > simi- 
larly indicating that M/r ~ |^ - ^FqI). The remaining physical 
possibility, po > Jo - 2/3o - 0, implies that ao = 0, which 
can lead to a non-zero finite limit for lim,^o Cr ■ ^' 
nite, then cr^ must not diverge, but if ^ is divergent (but not 
faster than logarithmic) as r — > 0, the exact limiting behavior 
of the corresponding cr^ should be inferred from the particular 
solution to the Jeans equations. 

By contrast, if M/r diverges (for which *!* — » -oo), then 
case (i) indicates that cr^ ~ M/ r — > oo whereas case (iii) re- 
quires cr^ > M/r and so cr^ — > oo (but ao = po = 0). In other 
words, cr^ necessarily diverges if M/r — > oo. Furthermore, 
cr^ must be divergent as fast as M/ r (note that if 6o > 0, then 
0-2 ~ M/r ~ 1^1, but 0-2 ~ M/r < ^ for 5o = where 6o is 
the negative logarithmic slope of M/ r or the potential) unless 
(>o - Jo - 2y6o = for which cr^ diverges faster than M/ r but 
not faster than I*!*!. 

For M ~ r (and *!* ~ In r), the result is basically that of 
equation dS); case (i) yielding the possibility of a finite lim- 
iting value of 0-2 whereas case (iii) is consistent with cr^ di- 
verging at most logarithmically or slower 

5.3. Infinite velocity dispersions? 

In the framework of classical Newtonian mechanics upon 
which the Jeans equations and the collisionless Boltzmann 
equation are ultimately based, the divergence of cr^. when 
M/r and the corresponding potential also diverge is in prin- 
ciple physically acceptable despite its mathematical quirk. 
However, it is clear that the arguments given in this paper 
eventually break down as o-^ approaches the speed of light. 
Moreover, in the corresponding halo, M/r should be diver- 
gent as r — > 0, and therefore there exists a radius below which 
GM(r)/c^ > r. Consequently, the central cusp, if it ever were 
present, must collapse to a singularity. In other words, one 
would expect that the formal infinity of the velocity dispersion 
can be always circumvented through the presence of a central 
black hole. The proper examination of physical behaviors of 
the tracers and the halo under these conditions would require 
consideration of relativistic physics, which is out of the scope 
of the current paper Of course, in reality, it is more likely 



that other various physical complexities in the system inter- 
vene to prevent the spherical Jeans equations to be applied 
uncritically all the way down to the center even before any 
relativistic eff'ects become important. 

6. PROJECTED QUANTITIES 

The direct measurement of radial and tangential velocity 
dispersions of stellar tracers is limited to nearby populations. 
More generally, the true observables are limited to the line- 
of-sight velocity dispersion - either the 'aperture-averaged' 
value or its profile for a subset. The implication of the the- 
orem on the behavior of the line-of-sight velocity dispersion 
is therefore of a great practical interest. However, we shall 
see that the integral transformation involved in the line-of- 
sight velocity dispersion weakens the theorem's practical con- 
straints. 

It is usually assumed that the observed line-of-sight veloc- 
ity dispersion follows the luminosity-weighted integration of 
the velocity dispersions along the line-of-sight direction. The 
latter ere is mathematically well-defined quantity such that 



(T 



2\ vcrlrdr 



(15) 



where 



Jr 



vrdr 



is the surface density of the tracers. If fieri is non-divergent, 
the leading term of Icr^ as 7? — > cannot be dominant over 
that of the surface density UK). Consequently, the leading 
term of o-^ in the central limit is largely dictated by the tracer 
density profile. 

In particular, if the density profiles of the tracers are approx- 
imated as power-law-like, we find the behavior o f the leadin 
terms for the surface density to be (see e.g. jAn & Zhaol200' 



V ^ Ar 



V ^ vo - Ar'' 



I ~ (j > 1) 

/~in/?-i (r=i) 
i~io-iiR^-^ (0<r<i) 

I ~ lo- hR^^'^ {Q<q<\) 
I ~ lo- IiR^lnR-^ (q^l) 
I~Io-hR^ (<?>!) 



where /q is the finite central surface density, and A and /i are 
some positive constants. Assuming v ~ r^'*' (y = if cored), 
0-2 ~ r" (a > 0), and fio > the corresponding behavior 
for equation ( fTSI l is similarly found to be 



lo-i 



\nR-^ 

Co + CiR"^^-y 
Co + CiR^ InR-^ 
Co + CiR^ 



(j > a + I) 

(y - a + I) 

(a - I <j < a + I) 

(y^a-l) 

(j < a - I) 



(16) 



with Co and Ci being some non-zero constants. Given the 
implication of the theorem for tracers with jo + 2fio in a non- 
divergent potential, i.e., o-r —> as r —> and so a > 0, 
we surmise that o-f ^ (~ /;™n(y-i,Q.)-) as 7? ^ if y > 1 
whereas it attains a finite limiting value (and typically increas- 
ing outward) if y < 1. If on the other hand (t^q is finite (for 
which Jo - 2/3o according to the theorem) or ySo = -oo (and 
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- cP^ non-zero), the leading term behavior of /cr^ is 
similar to that of / alone, and thus ct^q is finite. 

These essentially imply that the behavior of cr^ cannot in 
general be directly inferred from the leading term approxima- 
tion of cr^ alone, and that the strict constraint from the theo- 
rem is somewhat lost by going through the integral transfor- 
mation. Although one may deal with /cr^ instead of cr^ or can 
in principle invert the integral equation for cr^ to yield cr^ (as- 
suming some particular yS), this still indicates that inferring cr^ 
from cr^ involves analyzing higher-order behaviors of the lat- 
ter and thus requires high-precision measurements. Further- 
more, this is independent of the well-known degeneracy of cr^ 
and y6 in the inversion of cr^ in a sense that even though one 
possesses perfect a priori information on the uncertainties 
in the recovered cr^ are always amplified by inverting cr^. 

6. 1 . A central black hole 

The preceding discussion presumes the finite central poten- 
tial well, which is appropriate for the potential dominated by 
the halo that is cored or cusped not so steep as the SIS. If 
the potential however is dominated by the central point mass, 
the theorem indicates that cr^ ~ |^| ~ 1 /r and so /cr^ ~ R^^ 
for y > (i.e., cusped tracer populations) or /cr^ ~ InT? ' 
for y - Q (i.e., cored tracer populations). That is to say, the 
line-of-sight velocity dispersion of the population tracing the 
Keplerian potential is necessarily divergent with its logarith- 
mic slope, \d\vicP-Jd\nR\ being equal to min(l,y) where y 
is the 3-d density (negative) logarithmic slope of the same 
tracers unless the orbits of tracers are completely circularized 
toward the center. Nevertheless, the direct application of this 
inference to the observational results warrants caution since 
an assumption of 'infinite' resolution is implicit in the argu- 
ment. That is to say, the result is strictly relevant only if the 
observation can resolve the so-called sphere of influence of 
the central point mass. 

7. CONCLUSIONS 



I n this paper, we have established a general theorem - stated 
in § 15.21 - that makes precise the relationship between the cen- 
tral limiting values of the density and velocity dispersions of 
a stellar population, together with the potential. Our theorem 
gives all the mutually exclusive possibilities that can occur in 
a stellar system. We note that our theorem has straightforward 
applications to a number of astrophysical problems, including 
the kinematical modeling of the stellar populations in dwarf 
sphe roidal galaxi es and elliptical galaxies. 

In lEvans et al.l ('2009). we presented a simplified version of 
the theorem and argued that it is the consequence of the spher- 
ical symmetry assumptio n. However , after the appearance of 
the preprint version of .Evans et al.l (l2009l) . Scott Tremaine 
(priv. comm.) convinced us that the theorem is due to the 
'non-analytic' point at the center The spherical symmetry is 
of secondary importance and only indirectly responsible for 
the theorem by requiring a coordinate singularity at the center. 
The theorem in this respect might be understood as an incom- 
plete boundary condition imposed on the Jeans equations at 
the center resulting from the consideration of one-sided regu- 
larity. 

Applying to the real astrophysical problems, the true moral 
of our theorem is the urging of caution against interpreting re- 
sults based on simplifying assumptions. For instance, if one 
were to reconstruct the dark halo density from the observa- 
tions of the surface density and the line-of-sight velocity dis- 
persion profile of a tracer population, the seemingly benign 
assumptions of spherical symmetry and isotropy combined 
with a cored luminosity profile already severely restrict the 
possible halo density (it cannot be cusped!). Such idealized 
reconstructions are limited by the straitjacket imposed by the 
theorem, yet the restrictions may be non-existent in reality - 
not unlike assuming a spherical cow! 

The authors tha nk S. Tremaine for comments made on 
lEvans et^ (|2009l) . which led to the current work. The Dark 
Cosmology Centre is funded by the Danish National Research 
Foundation. 
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